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ABSTRACT 


In the classical gray atmosphere, the phenomenological 
coefficient of absorption that enters into Kirchhoff s law is 
assumed to be independent of the frequency ( k v = const = *). 
Strictly considered, this coefficient is the product obtained 
by multiplying the Einstein coefficient of true absorption by 
a frequency- dependent correction for induced emission 
(Rosseland factor). The alternative here examined is to set 
the Einstein coefficient constant, incorporate the effects of 
induced emission into the transfer formalism, and deter- 
mine the march of the local thermodynamic equilibrium 
(LTE) source function by numerical iteration. This second 
kind of gray atmosphere differs markedly from the classi- 
cal one in respect to many physical characteristics and 
represents, in conventional terminology, the simplest non- 
gray problem that is physically realistic. 
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RADIATIVE EQUILIBRIUM IN AN ATMOSPHERE WITH 
CONSTANT EINSTEIN ABSORPTION COEFFICIENT 

by 

Rupert Wildt* 
and 

Sandra Schwartz* 

Goddard Institute for Space Studies 


INTRODUCTION 


The theory of radiative transfer through the gray atmosphere antedates the concept of induced 
emission and the ensuing distinction (Reference 1) between two kinds of phenomenological coeffi- 
cients of absorption [dimension (cm' 1 )]. The first of these is properly called the Kirchhoff coef- 
ficient, because it appears in Kirchhoff 1 s law, and the assumption that it be gray ( k v = const = k) 
is fundamental to the classical theory of k. Schwarzschild (Reference 2) and its elaboration by Milne, 
Hopf, and others. The second one (Reference 3), which measures the attenuation by true absorption 
but excludes the effect of induced emission (negative absorption), will here be denoted by and 
called the Einstein coefficient. Owing to the exclusion noted, the inequality k v < a v holds generally. 
The equation 


K 


v 



(1) 


while true for thermodynamic equilibrium (Reference 1), is unnecessarily restrictive (Reference 
4). A relation of unrestricted validity, 


K 


V 



( 2 ) 


where (erg/cm 3 -sec-ster) is the phenomenological coefficient of spontaneous (isotropic) emis- 

sion, follows from the most general, time- dependent form of the equation of energy transfer in 
Cartesian coordinates, 

i c 2 i v 

c <5T + Mi -giq + m 2 + m 3 = ~ a u I u + e u + -^r e v - (3) 


*Yale University Observatory 



III 


1 



in which the first term on the right-hand side represents the rate of true absorption, the second 
term the rate of spontaneous emission, and the third term the rate of induced emission, which must 
occur precisely in the direction of the incident pencil and at a rate proportional to its intensity, ac- 
cording to Einstein’s analysis of the momentum transfer. By rearranging the terms on the right- 
hand side of Equation 3, the conventional form of the energy transfer equation is recovered, 

i*K d K 

c + Mi <9 Xl + M 2 (9 x 2 + M 3 dx 3 

and the relation (2) is thus seen to apply even to non-isotropic and transient radiation fields. 

A more transparent form can be obtained by introducing the familiar definition of the source 
function, 



and writing 



(5) 


( 6 ) 


Instead of using i* 9 the spectrum of the source function can be characterized by a reciprocal mono- 
chromatic temperature, 


1 




(?) 


which will generally be frequency- dependent. An equivalent version of Equation 6, then, is 


K 


V 



( 8 ) 


The great interest attaching to Equations 6 and 8 derives from the fact that they rest on the strictly 
formal nature of the definitions 5 and 7 which gives them validity for any spectral distribution of 
the source function. 

In the special case of LTE, Kirchhoff T s law holds, by hypothesis, so the right-hand side of 
Equation 5 becomes the Planck function: 


I 


* 

V 



h.v 

exp-^ 




- f 

Jo 


I * dv , 


(9) 


2 


where <? is the Stefan- Boltzmann constant. In this case, both sides of Equation 7 are independent 
of the frequency, 


1 




( 10 ) 


and Equation 8 reduces to Rosseland’s relation, (Equation 1). If, in addition to the postulate of 
LTE, it is assumed that the Kirchhoff absorption coefficient is independent of the frequency (classi- 
cal gray atmosphere with k v = const = *), then it follows necessarily from the relation 


K 



(ID 


that the concealed mechanism of induced emission implies 


lim - 00 (12) 


and that the rate at which the Einstein coefficient tends to this limit varies from level to level in 
the atmosphere, depending on the local value of i*. As in the classical model of a semi-infinite at- 
mosphere, I* itself tends to infinity with increasing depth; similarly, at any finite frequency, the 
Einstein coefficient, a v , tends to infinity. The presumed simplicity of the classical concept of 
gray matter turns out to be specious when due attention is paid to the phenomenon of induced emis- 
sion. Conceding that the classical gray atmosphere is a rather artificial model detracts in no way 
from the value of the study lavished on it. It does suggest, however, examination of a second kind 
of gray atmosphere, distinguished by postulating the Einstein coefficient to be independent of fre- 
quency; in other words, let it be assumed that 


- constant - a 


and 


K 


V 



or 


K 


V 



(13) 


(14) 


(15) 


Radiative transfer through such an "Einstein- gray" atmosphere poses, perhaps, the simplest non- 
gray problem (in conventional terminology) that is physically realistic. 
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The well-known formulary of non- gray radiative transfer in a plane parallel atmosphere is 
readily adapted for use in the Einstein- gray case by introducing as independent variable the 
neutral optical depth defined by 


r 


r 

■'x 




(16) 


Moreover, by Equation 14, 



(17) 


and the optical depth that measures the monochromatic attenuation of a pencil escaping from level 
t below the surface is 


(r) 


- f 


1 + 


(t) 


2h v 3 


dt 


(18) 


Then the condition of strict radiative equilibrium at the neutral level reads as follows: 

C(t)E, [t„(t)-T„ (T)} dt 


f 

c’CCrj] 

“J 

L, , * r 

c 2 C(t) 

Jo 

1 + , 

2hv 3 J 

1 


1 + , , 
2hv 3 


-if 


1 + 


C 2 I * (t) 


2hv 3 


I v* (t) E i [ r v ( T ) " t v (*)] dt f = o - 


L&Lll 

4 dr » U9) 


where F(r) is the net radiant flux at the level r. This is the analogue of Milne T s integral equation for 
the classical gray atmosphere. The use of the symbols I* or T* in Equations 14 through 19 stresses 
the fact that they hold in the absence of LTE. 


THE EINSTEIN-GRAY ATMOSPHERE IN LTE 

When the Planck function (Equation 9) for I* is substituted in Equation 19, the first term in- 
side the braces (multiplied by the preceding Rosseland factor) can be made integrable over the 
entire spectrum, so that the integral is directly proportional to the total source function I* (r). 
This method provides no such simplification for the second and third terms, because the mono- 
chromatic optical depths appearing as arguments of the exponential integrals, E 1? preclude separa- 
tion of the integrations with respect to v and t. Since this LTE problem is non- gray in the con- 
ventional (Kirchhoff) sense, there is small hope of finding an exact analytical solution. Therefore, 
it has been thought worthwhile to try an iteration scheme in order to obtain a numerical solution. 
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After reverting to the familiar symbol B : , (r) for the LTE source function, and writing 


\ (t) 



B(r) 


f 

J o 


B„ (r) dv , 


T(r) = [vB(r)/of 


( 20 ) 


the problem to be solved is to find the total source function B(t) satisfying the integral equation ex- 
pressing strict radiative equilibrium, or constancy of the net flux ttF of the total radiant energy, in 
a gray atmosphere of the second kind, namely, 



in which B v (r) is given by Equation 20, and j v (r) is the mean monochromatic intensity. The latter 
is related to B v (r) by Hopfs A- operator. Since the Rosseland factor appears in the monochromatic 
optical depths entering the A- operator, its explicit form in this case is 


M b , 



j i/ ( t ) • 


with 



( 22 ) 


(23) 


and 



The condition of flux constancy, 


(24) 


F(t) 



(r) dv 


constant , 


0 < r < oo , 


(25) 


5 




is equivalent to the integral Equation 21. 

This step completes the statement of the problem. Its numerical solution will be carried out 
by adapting to the problem at hand an iteration scheme due to Lecar (Reference 5). 


THE LINEAR APPROXIMATION TO THE SOURCE FUNCTION: B (0) H 


The first step is to find a linear approximation to the source function that will yield the asymp- 
totic value at large optical depths of the prescribed total net flux. 

As Lecar has proved, by expanding (t J in a Taylor series about the point t v - r v , 


, . 4 , * B , K) 

limF^r) = 


4 dB v dT 3t 

~ST WST- v - 


(27) 


provided that 


lim = 0 , 

r -00 St n 


n > 3 . 


(28) 


According to Equation 24 




(- 


hv \ - 1 
kT(r)j 


(29) 


Therefore 


dT 

r | 

- 

/ M" 

dr 

<9T 

1 - exp 

(” kT j 


Jo 

_ 



- F 


(30) 


where F is the prescribed value of the total radiative flux. 
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With k r defined by 


1 dB 
k r dT “ 



Equation 30 becomes 


± , . _l_dBdT _ 4_ , _1_ dB 

3 k dT dr “ 3 k dr “ F • 

co r t — > 00 r 


(31) 


(32) 


Substituting in Equation 31 the following expressions, 


277 5 k 4 
15c 2 h 3 ’ 


dB 

dT 


3 

15c 2 h 3 


and 


d JK 2h 2 v 4 

dT c 2 kT 2 


exp 


(w) [«P (w) - .] 


shows that 


K = 4tt 4 k 3 t* I 1/4 exp (w) [ exp (w) - l] 


dz/ 


In terms of the variable x, defined by 


(33) 

(34) 


(35) 


(36) 


hu 

x - kT 


and 


dx 



(37) 


(38) 


1 

k 

r 


15 f” x 4 e 2x 

4tt 4 J 0 (e«-l) 3 


dx . 


(39) 
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Integration by parts shows that 


f 

J o 


(e x - l) 3 


dx - 6 


r x 2 J r x 3 

Jo e * ~ 1 X + 2 J 0 e x - 1 


dx 


(40) 


The integrals on the right side of Equation 40 can be expressed in terms of the zeta function £(s), 
where 


Therefore 


and 


£(») 


y*C0 00 

L*- 


(41) 


00 00 




(42) 


, n = 1 n= 1 


k = 0.94758593 


(43) 


Let the desired linear approximation to B(r) be denoted by B< 0) (r) . Then, in order to obtain 
the correct value of the total flux at large optical depths, B <0) (r) must satisfy Equation 32 every- 
where. Therefore for all r 


dB (0 > (r) _ 3_ 

dr 4 


k F 


(44) 


and B<°>, (r) has the form: 


B<°> (r) = -fk r Fr + B<°>(0) . 


(45) 


The constant B<°> (0) is arbitrarily fixed by requiring the function B (0) (r) to satisfy the equation 


2 


Jo 


B<°> (t) E 2 (t) dt 


F . 


(46) 
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If the problem under discussion were the classical gray atmosphere, Equation 46 would stipulate 
that B (0) (r) yield the prescribed net flux at the surface. This physical interpretation of Equation 
46 does not hold for the gray atmosphere of the second kind, as is shown by Equation 26. The re- 
sult of substituting the expression for B (0 > (r) from Equation 45 in Equation 46 is 

B< 0) (0) = \ (2-k r ) F - 0. 526207F . (47) 


Therefore the linear approximation to the source function is 

B<°> (r) = 0.710689F[r + 0.740418] . (48) 

THE SECOND APPROXIMATION TO THE SOURCE FUNCTION: B ( 1 ) (r) 

The second approximation to the source function, B (1) (r) , is obtained by making two correc- 
tions to the linear approximation B (0) (r) . The two corrections are denoted as h t B (0) (r) and 
8 2 B (0) (r ) , and 


B (1) (r) - B (0) (r) + Sj B (0) (r) + S 2 B (0 ) (t) 


(49) 


The first correction, B c °> (r), is defined in terms of the deviation from F of the integrated 
flux associated with B (0) (r) . Finding h 1 B (0) (r) requires an explicit expression for the radiative 
flux that would exist in the gray atmosphere of the second kind if B <0) (r) were the source 
function. 


The expression for the monochromatic optical depth from Equation 23 will now be used to 
transform Equation 26 into an expression for F v (r) which no longer contains the variable r v , i.e., 


F u (r) - 2 <1 - exp 


{ | - ex P [" kT(i)]} (t)E a jl - exp [- k T(^] 


d <? dt 


“ 2 J o j 1 " exp [" kT(t)]} (t) E 2 j 1 - cxp [" 

On replacing B v (t) by 2h^ 3 /c 2 jexp [hv/kT(t)] - i| \ it follows that 


d<f ] dt „ (50) 


4hv 3 

f *( t ) = -pr 


l " p [■ «<•>] e > (I {■ ■ [■ 


i\ 


d«f] dt 


4hv 3 

~ 2 


kT(f )J J 

I exp [■ (l { ~ exp [- kfe]} d<f ) 


dt . (51) 


9 


To simplify this expression, write hv/kT(t) in the form 


T 

hv hv 1 e 

kT(t) ” kT e T(t) * 


( 52 ) 


where 


T 



(53) 


is the effective temperature of the atmosphere. Then define the variable y and the function z ( 7 ) by 


hu 

y - kt 


(54) 


and 


Z ( T ) " T( 


t) [B(r)J 


(55) 


so that 


hv 


kT( t ) 


- yz(t) . 


(56) 


In terms of these variables 


( T ) 


^F^y 3 j e ' yZ<t)E 2 |j r [l-e-^>] 


d*>dt 


- F w; y 3 J e ' yz(,) E 2 |j t D - e ' yz(f) ] dt • 


(57) 


With f (r) defined by 


f y ( T ) 


™ y 3 | [l - e-^( 0 ] 


dt 


~ y 3 | e- yz <‘>E 2 [l - e" yz ^>] dt , 


(58) 
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(59) 


K CO 



The integrated flux at the level r is 


F(t) 


F 



0)dy . 


(60) 


The integrated flux that corresponds to the source function B<°> (r) will be referred to as 
F (0) (r) . It can be found by replacing B (r) by B<°> (r) in the expression for f y (r) given by Equa- 
tion 58. Accordingly, set 


z< 0) ( T ) 



(61) 


and 





• / T 


e -yz ( 0 ) (t> g 




dt 


r, 

w Jo 


•yz ( 0 ) (t) j. 



e - yz (0) (£) 


^ dt . (62) 


Then 


F (0) (r) = F f y <°> (r)dy . (63) 

•'o 

The extent to which F (0) (r) deviates from the prescribed net flux F is measured by the func- 
tion F (0) (r), where 


8F<°) (r) = F - F (0) (r) 


•^o 


f y (0) (r) 


dy 


(64) 


The correction h x B<°> (r) is defined in terms of SF (0) (r) by an equation similar to Equation 44, 
which defined B<°> (r) in terms of F, i.e., 

= |k r SF<°>(r) - (65) 


11 


The solution to this equation is 


S j B (0 ^ (r) = |-k r j SF<°> (t) dt + 5jB (0 > (0) . (66) 

Jo 

The constant b (0) (0) is chosen in terms of SF (0) (0) in the same manner that B (0) (0) is defined 
in terms of F in Equation 47; i.e., 


8 1 B<°>(0) = \ (2 -k r ) SF<°> (0) . 


(67) 


The first correction to B (0) (r) is then 


SjB(°)(t) = -fk r f SF (0 > (t) dt + 8F (0 > (0) (2-k r ) 

Jo 


( 68 ) 


or, with k r = 0.94758593, 


f 


(r) = 0.710689 SF<°> (t ) dt + 0. 526207 8F<°> (0) 


(69) 


The second correction to B (0) (r) , S 2 B (0) (r) is defined in terms of the derivative of the in- 
tegrated flux associated with B (0) (r) . This derivative will be referred to as dF(0)/dr. It can be 
expressed in terms of B (0) (r) by rewriting Equation 21 so that all the terms are explicit functions 
of B<°> (r). 

The monochromatic mean intensity J v (r) according to Equation 22 is 


h [ r v ( T )] 



(70) 


Since 


( r ) - 



dt 
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and 


r t 2hz/ 3 , 

K Mo] = -rr i ex p 




hv 3 

Jv (^) = 


j exp k^y] E, ( £ |l - exp ^ kT h (^) 



dt 


(71) 


In terms of the variables y and z of Equations 52 and 53, 


7.5 h 3 

J v (r) = ^r F vr y 


J e-y^E, j [l - e- yz ^>] d£ 


dt 


(72) 


Let a quantity j (t) be defined by 


j v ( T ) 


r 00 

^■1 


Ej "S II [l " e' yz ^ 5 ] d£ 


dt 


(73) 


Then 

Jv( T > = F W7 jy( r > • 

The function (r) can be written as 

b„(t) = kY; y 3 [e y2(T) - I]" 1 • 

With b y (r) defined by 

b,(r) = ^ 1 e,!(,> - A"* ’ 

B„(t) = FkY7b y (r) . 


(74) 


(75) 


(76) 


(77) 
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When these expressions for J„ and B v are substituted into Equation 21, it becomes: 


= F [l - [j y (r) - b y (r)] dy = 0. (78) 

•'O 

The function (1/4) (dF/dr) <°> can be found by replacing B(r) with B<°> (r) in all the terms on 
the right side of Equation 78. In order to do so, set 

j y <°>(r) = fy 3 _[ e-^ (0) <‘)E 1 ||| [l - e-^ <0) tf)] df|| dt (79) 

and 

b (0) (T) = y 3 [ey‘ (0) < T > - lT‘ . (80) 

y 77^ L J 

Then 

\ (^) <0> = F j [i-e-^ <0) <->] [j y <°> (r) - b y <°) (r)] dy . (81) 

This function is not identically equal to zero because B (0) (r) is only an approximation to the source 
function. 

When B (0) (r) is corrected by S 2 b (0) (t), the monochromatic Planck function associated with 
B<°) ( T ) is also corrected. This monochromatic Planck function is B/ 0) (r), where 


B/ 0) (r) 


kT, 


b y <°> (r) 


(82) 


The correction to B/ 0) (r) associated with S 2 b (0) (t) will be referred to as S 2 B/°> (r) . It is de- 
fined by the following equation: 



exp 


hv 

kT' 0) (r) 


S 2 B/ 0) (t) dy 



(83) 


where 


T<°> (r) 


ttB< 0) (r )j 4 


(84) 
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Equation 83 can be transformed into an explicit equation for S 2 B (0) (r) by considering the changes 
in the temperature distribution within the atmosphere associated with the corrections S 2 B„«» (t). 

If S 2 T (0) (r) represents the temperature corrections corresponding to S 2 b (0) (t), then 


S 2 B v < 0 )(t) = B„[t<°>(t) + S a T(°)(r)] - B„ [T<°>(t)] . 


(85) 


Expanding B v [t ( 0) (r) + S 2 T (0 > (r)J in a Taylor series of powers of S 2 T (0) (r) gives the following 
expression for S 2 B J 0) (r): 


00 

= ££ [«,*•> (r)]' 


d" B„ (T)~ | 
<?T n 


J T=T„ 


( 86 ) 


If all powers of S 2 T <0) (r) except the first are neglected, this expression becomes: 


S 2 b„<°>(t) = s t<°>( 


<?B„ (T)l 

dT J. 


(87) 


The correction function S 2 b ( 0) (r) can be expressed as: 

S 2 B<°>(t) = b[t<°> (r) + S 2 T(°) (r)] - b[t<°> (r)] . (88) 

If b[t<°> (r) + s 2 T<°> (r)] is expanded in a Taylor series in powers of s 2 T (0) (r), Equation 88 
becomes 


S 2 B(°> (t) 


ZZ^( S 2 T(0) (^)) n 

n = 1 


d n B(T) 
” dT n 


T= T 


(0) * 


(89) 


By keeping only the first terms of this series: 

S 2 B<°>(t) = 

It follows from Equations 87 and 90 that 

S 2 B y <°> (r) 
8 2 B<°> (t) 


: T( ° , (-)#] t=t( o ) - 


dTJ T=T (0) 

dB~ | 

dT J T =T <0) 


(90) 


( 91 ) 
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and therefore 


S 2 B(°)(t) = 



S 2 B (0 > (r) 


(92) 


Substituting this expression for S 2 Bj- 0) (t) in Equation 83 gives the following equation for 


S 2 B (0) (t): 


^0 


1 - exp 


hv 


kT< 0 > (r) 


fK 1 

dTj t=t (°) 1 

S 2 B (0 > (t)\ 5b _ I dv 

^(9T_ T=T (0)y 


1 /dF\ (0) 
~ 4 \dr ) 


(93) 


With k defined by 


, dB 

ft 

PQ 

8 

■ 

( hzA 

k G dT ” 

) (9T 

J o 

1 - exp 

tw)J 


dv 


(94) 


Equation 93 becomes 


... 1 /dFV°> 

8 2 B<) < t > = 4k^dF] 


(95) 


From Equations 33, 34, and 35 it is found that 


15 h 5 
4w 4 k 5 T 5 


f 1,4 ex p(^f) ~ 1 

•'o L N 


dv . 


(96) 


If the quantity hv/kT is replaced by the variable x, the expression for k G becomes 


15 
4 tt 4 


r x 4 

Jo ^ - 1 


dx . 


(97) 


This integral can be written in terms of the zeta function £(s) of Equation 41: 


CD 



k= 1 


0.958057 . 


(98) 
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The second correction to B (0) (r) is therefore 


s 2 B<°) (r) 


0. 260945 



(99) 


The second approximation to the source function is found by adding the correction in Equations 
69 and 99 to B<°> (r); i.e., 


B (1) (r) = B<°) (r) + 0. 710689 SF <0 > (t) dt + 0. 526207 SF (0 > (0) 

Jo 


/dF) (0) 

+ 0. 260945 


( 100 ) 


THE METHOD FOR FINDING B (i+1 >H FROM B (i) (r) FOR i >:0 

The procedure used to obtain B (1> (r) from B <0 > (r) is the basis of the iteration scheme that 
gives successively higher approximations to B(r). Let us assume that B (i) (r) is known. Then for 
i >0, 


B < i+1 )( T ) = B (i 5 (t) + SB< i > (r) 


( 101 ) 


where 


f 

Jo 


SB (i) (r) = 0.710689 SF^ >( t ) dt + 0 . 526207 SF^ 1 > (0) 


+ 0.260945 


/dF\ (i) 
\dr / 


SF (i) (r) = F - F (i) (r) , 


F (i) (r) = F f<‘>(T)dy 


r 

Jo 


( 102 ) 

(103) 

(104) 
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and 

b y (i) (r) = ^y 3 [e yz( ° (T) - l] (109) 

NUMERICAL METHODS USED IN CALCULATING THE SET OF 
FUNCTIONS B (i) (r), i = 0, 1, • • • , 8 

The functions B (i) (r) and §B (i > (r) depend on the value of the prescribed flux F. This is not 
true of B (1) (t)/F and SB (1) (r)/F. Therefore, B (1) (r)/F and SB (1) (t)/f can be computed without 
specifying F. 

Using the IBM 7094 and 360-75, the functions B (i) (r)/ F and 5B (i) (r)/ F, i = 0, 1, • • • 8, have 
been evaluated at the following 652 values of r: 

Tj = 0 , r 2 - 10“ 4 , * * ■ , 7" n ~ io- 4+ ' 01 ( n " 2 ) , • • • , r 602 = 100 , • • * , 

T 65 2 ~ 316 . (110) 

The following additional set of equations was used to compute the monochromatic quantities 
f y (i) ( r j) and j y (i) ( T j) at j = 0, • * • , 652, and i = 0, 1, • • • , 8: 

kT 

f y (i) ki) = (ill) 
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F y < 0 fr) = 2 f (i) , , B y (i> (*y) E * fry ' -y^fi)] dt y 

^ T y 


- 2 


r T y ( ‘ ) ( r i) 

J o B y (i) ( t y) E 2 [% (i) (r,) ~ t y ] dt y ; (112) 


(i) 


(T ‘> = !.'( 


exp <- y 


B (i 


i) (T) | fl dt y 


kT 


j y <0 ( T i) = “if J y (i) h)/F ; 


(113) 


(114) 


J "’ ( T i) = T J o »,“> (*,) E, [r„<‘> (r,) -.Jdt, 


-00 

tJ M E > [', - -,<■> (-,)] »■, ■ (115) 

*% ( T i) 


kT 


b y (1) ( T i) = “f B y (i> ( T y)/ F ! 


B y (i) ( T i) = ;t f wt, y 3 l ex P 



(116) 


(117) 


The exponential integrals E n (x) were approximated by a sum of 10 terms by means of the 
Gauss- Legendre quadrature formula. By definition 


£ n (x) - jf K n 2 exp (- d\ . 


(118) 


In terms of the variable 


f = 2\ - 1 


E n< X ) " | j (£ + X > n ' 2 eX P (" fTf) 


df 


(119) 

( 120 ) 


19 



I 


II 


II I 1*1 ■ 


The Gauss- Legendre formula states that 

E nO) = ^7 [ (£ + l)"~ 2 exp jE^riJ dg = ^ A i(^i +1 )"' 2exp (r FTTjj ■ (121) 

In Equation 121, the 10 constants are the roots of the equation 

Pio (x) = 0 , (122) 

where P 10 (x) is the 10th order Legendre polynomial, and the 10 constants A ; are 

A i = [Pio ( x i )]" 2 (irfr) . i = 1 . • • • 10 . ( 123 ) 

In the process of choosing values of y, the y-axis was divided into intervals of length 5. Within 
each interval, 10 values were chosen, corresponding to the 10 zeroes of the 10th order Legendre 
polynomial, normalized for that interval. For i = 0, * • • , 8, the largest value of y, , was 
chosen to satisfy the equation: 


m T ax l j y (i> ( T i) “ b y ( ° ( T i)l < 10 " 8 


(124) 


where 


y > , j = 0 , * * • , 652 . 


The monochromatic quantities f y (i) (r.) and j y (i) (r.) were computed from Equations 111 
through 117, with t as the depth variable. The procedure followed was similar to the one used by 
Lecar (Reference 5), except that the functions E n (x), n = 1, 2, were evaluated by means of Equa- 
tion 121. It was assumed that b (r) is a linear function of r for r > r 6S2 and that 


B y (i) (r) 



(125) 


where 


T > T 


65 2 


The integrations in the variable y, necessary for the evaluation of SF (i) (t.) and dF (1) /dr [r . ) f 
were carried out by assuming 

f y (i) ( T i) = jy (i) ( T i) " b y (i) ( T j) = °’ ( 126 ) 
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where 


y > 4(*> . j = 0, * • • , 652 


and by using the Gauss- Legendre integration method for 0 <y <>£(*) . 
The integrals 


0.710689 



$F (i) (t) dt 


were evaluated by means of the trapezoidal rule. 


THE CONVERGENCE OF THE SET OF FUNCTIONS B (i) (r)/F, i=0, 1, ■ ■ ■ , 8, 

TO THE SOURCE FUNCTION B(r)/F 

The ideal outcome of the use of this iteration scheme would be the convergence of the set of 
functions B (1) (r)/F to a function B (k > (t)/f for which 


SF °° ( T n)/ F = 0 • 

(127) 

/ dF \ k . . , 
idfj K)/ f = o , 

(128) 


8B<k) ( r n)/ F = °- n = 1, 2, • - ■ , 652 . (129) 

The function B (k) (r)/F would then be exactly equal to the source function B(r)/F at r = r l9 • • • , t 652 . 

An indication of the convergence of the computed sequence of functions B ( i} ('O/F, i = 0 , 1, ■ • • , 8, 
is provided by Table 1. Values are listed for r n < 100. 

Table 2 gives a more complete indication of the deviation of F (8) (t)/f from 1 and the deviation 
of dF/dr (8) (r)/F from 0. The numbers in Table 2 have been obtained by interpolation from the val- 
ues of F< 8 > (tJ/f and dF/dr( 8 > (tJ/f, n = 1, • • • , 602. 

Tables 1 and 2 indicate that the function B (8) (r)/F comes close to satisfying Equations 127 to 
129. Therefore the numerical values computed for B (8) (r)/F can be used as a sufficient represen- 
tation of the source function B(t)/f. 
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Table 1 


Quantities Illustrating the Convergence of the Sequence of Functions B ( 1 5 (r)/ F, i = 0, 1, • • • , 8. 


i 

The maximum 
value of 

8B (i) (r )/F 

The value of r at 
which SB< i > (r)/F 
has a maximum 

The maximum 
value of 

SF^> (t)/F 

The value of r at 
which 5F< i > (r )/F 
has a maximum 

The maximum 
value of 

l/F (dF/dr) (i) (r) 

The value of r at 
which l/F (dF/dr) (i) (r) 
has a maximum 

0 

-0.0783748 

0.000000 

-0.020569 

0.000000 

0.300350 

0.000000 

1 

-0.0050028 

0.042658 

+0.003884 

0.501187 

0.017953 

0.042658 

2 

-0.0007000 

0.083176 

+0.000711 

0.512861 

0.002573 

0.085114 

3 

+0.0001200 

1.122018 

+0.000134 

0.524807 

0.000427 

0.117490 

4 

+0.0000247 

1.122018 

+0.000030 

0.524807 

0.000325 

99.999985 

5 

+0.0000051 

1.071519 

-0.000005 

2.041738 

0.000325 

99.999985 

6 

+0.0000011 

1.096478 

-0.000003 

1.995262 

0.000325 

99.999985 

7 

+0.0000002 

1.148154 

-0.000003 

2.137962 

0,000325 

99.999985 

8 

+0.0000000 


-0.000003 

2.187761 

0.000325 

99.999985 



VALUES OF BM/F AND THE REMAINDER 
FUNCTION rW 

Conventionally, the solution of Milne f s in- 
tegral equation for the classical gray atmos- 
phere is represented as the sum of a linear 
term and a bounded remainder function, viz. 

B(r) = f f[t + q(r)] . (130) 

The variable r in this solution is defined in 
terms of the absorption coefficient given by 
Equation 11. But the same symbol r is used 
throughout this paper for another variable de- 
fined by Equations 13 and 16. To prevent con- 
fusion in what follows, the conventional notation 
of Equation 130 will here be changed to 

i b( i t) = “4 F [ T r + q( r r)J . (131) 


Table 2 

The Deviation of 1/FF (8) (r) from 1 and the 
Deviation of 1/F (dF/dr) (r) from 0. 


T 

1/FSF <8 > (t) 

l/F (dF/dr) < 8 > (r) 

0.00 

+0.0000005 

-0.000000 

0.01 

+0.0000005 

-0.000001 

0.05 

+0.0000005 

-0.000001 

0.10 

+0.0000005 

-0.000001 

0.50 

-0.0000005 

+0.000001 

1.00 

-0.0000017 

+0.000000 

1.50 

-0.0000024 

+0.000003 

2.00 

-0.0000026 

+0.000007 

5.00 

-0.0000021 

+0.000026 

10.00 

-0.0000018 

+0.000052 

25.00 

-0.0000015 

+0.000121 

50.00 

-0.0000011 

+0.000212 

75.00 

-0.0000008 

+0.000278 

100.00 

-0.0000006 

+0.000325 


This notation emphasizes the physical distinction between the independent variables 1 r and r, 
which are not strictly comparable. 

It is convenient to represent the solution for the source function in a gray atmosphere of the 
second kind in a form similar to Equation 131, writing 


B(r ) * 0.710689 F[r + r(r)] . (132) 

Some numerical values ofB(r) and of the remainder function r(r) are listed in Table 3. Evidently 
the source function B (r) increases monotonically for r > 0. The remainder function r(r) is bounded, 
is always positive, and attains a maximum near r = 1.5. 

If instead of 0.710689 another constant, A, had been chosen on the right-hand side of Equation 
127, i.e., 


then 


B(r) = AF[r + s(r)] , A / 0.710689 


A). 710689 0.710689 , x 

*00 = { A V r + ~A r ( T > 


(133) 


( 134 ) 
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and s(r ) would be unbounded. The preferred 
constant 0.710689 is identical with the constant 
appearing in the linear approximation to the 
correct source function (Equation 48). The 
reason for this identity is not obvious and bears 
further investigation. 

In order to estimate the accuracy of the 
solution found for B(r) and r (t) , the functions 
X B ( x r) and q ( x t) were computed in a similar 
manner, using Lecar's iteration scheme and the 
numerical approximations of the earlier section 
of this paper, ’’Numerical Methods." The values 
found for q(0) and q( x ) were correct to six 
decimal places. The maximum deviation from 
F in the net flux associated with this computa- 
tion occurred near r = 0.1 and was equal to 
10 " 6 . The maximum deviation from 0 in the 
derivative of the net flux occurred near r = 1 
and was equal to 4 * 10"° . A comparison of the 
magnitude of these errors with the errors dis- 
played in Table 2 suggests that the values of 
B(t) and r(r) listed in Table 3 are probably 
correct to five decimal places. 

THE SURFACE TEMPERATURE 

An exact result of the theory of the classi- 
cal gray atmosphere is the Hopf-Bronstein 
relation: 

Iq(0) = 7f ' dss) 

Hence, by Equation 131, 

r B(0) = "^-F . (136) 

Then from Equations 20 and 53 it follows that 
the boundary temperature is 
1 

I T(0) = T e = 0. 8112 T e . (137) 


Table 3 


Values of the Source Function B(r)/F and the 
Remainder Function r(r) . 


T 

r(r) 

B(r )/F = 0.710689 [r + r(r)] 

0.00 

0.62418 

0.443597 

0.01 

0.63486 

0.458293 

0.02 

0.64253 

0.470853 

0.03 

0.64855 

0.482240 

0.04 

0.65360 

0.492937 

0.05 

0.65803 

0.503191 

0.06 

0.66202 

0.513133 

0.07 

0.66568 

0.522837 

0.08 

0.66906 

0.532348 

0.09 

0.67221 

0.541695 

0.10 

0.67516 

0.550898 

0.20 

0.69700 

0.637489 

0.30 

0.71065 

0.718260 

0.40 

0.71999 

0.795964 

0.50 

0.72667 

0.871779 

0.60 

0.73155 

0.946318 

0.70 

0.73516 

1.019951 

0.80 

0.73783 

1.092918 

0.90 

0.73979 

1.165383 

1.00 

0.74122 

1.237466 

1.10 

0.74222 

1.309249 

1.20 

0.74290 

1.380799 

1.30 

0.74332 

1.452164 

1.40 

0.74353 

1.523385 

1.50 

0.74358 

1.594490 

1.60 

0.74350 

1.665502 

1.70 

0.74332 

1.736442 

1.80 

0.74305 

1.807322 

1.90 

0.74273 

1.878158 

2.00 

0.74235 

1.948957 

5.00 

0.72940 

4.071821 

10.00 

0.72060 

7.619014 

25.00 

0.71459 

18.275088 

50.00 

0.71237 

36.040746 

75.00 

0.71152 

53.807382 

100.00 

0.71107 

71.574293 
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The corresponding relations for the gray atmosphere of the second kind are 


r(0) = 0. 62418 (cf. Table 3) 


(138) 


and 


B( 0) = 0.710689 r(0) F = 0.44360F, 


( 139 ) 


from which it follows that 


T( 0 ) = 0. 8 161 T e . 


( 140 ) 


The continued inequality 


: T(0) < T(0) < T c 


( 141 ) 


calls for some comment. 

The effective temperature can be regarded as that of an isothermal blackbody, of infinite 
optical thickness, whose surface emits the same flux as emerges from the two atmospheres here 
compared. The inequalities stated suggest that the temperature stratification in a gray atmos- 
phere of the second kind is closer to isothermy than is that in a classical gray atmosphere, or, in 
other words, that the average temperature gradient throughout the layers substantially contributing 
to the escaping radiation is smaller than in the classical case. A precise comparison of the two 
temperature gradients is impossible because the optical depths r and l r differ with respect to their 
physical definition. Consequently, the scales of r and l r cannot unambiguously be related to a 
common geometric scale. In these circumstances a somewhat arbitrary, though plausible, as- 
sumption has to be made. Let t* and x r* be the optical depths at which the respective source func- 
tions equal the flux constant. Then 

r* = 0.6728, B (t* ) = F, (142) 


and 


x r* = 0.6454, i b( i t*) = F. 


( 143 ) 


It is assumed that these levels are approximately equivalent physically, in the sense that the 
average temperature gradients between the starred levels and the surfaces lend themselves to a 
meaningful comparison. Now these gradients are 



( 144 ) 
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and 



4.4882 . 


(145) 


The gradient of Equation 144 is indeed smaller than the classical gradient (Equation 145). 


THE MONOCHROMATIC FLUXES: F v [ t ) 

The frequency distribution of the radiant fluxes at any level in the atmosphere is given by the 
function ttF v (t) dv [erg/cm 2 -sec] appearing in Equation 26. According to Equation 59, the normalize* 
monochromatic net flux 


f y (?) 


F v (r) kT^ 
F h 


(146) 


is a function of r only. Some values of f y (r) axe given in Table 4. They were found by interpola- 
tion from the more extensive set used in the computation of B(t)/F. 

The frequency distribution of the flux changes progressively as the radiation flows towards the 
surface. The monochromatic flux at high frequencies decreases, while the monochromatic flux at 
low frequencies increases. The frequency at which the monochromatic flux is a maximum shifts 
towards the red end of the spectrum, and the maximum increases as the radiation approaches the 
surface. This pattern of redistribution, or flux reddening, is illustrated in Figure 1, which con- 
tains graphs of the function f y (r), evaluated at r = 4, 3, 2, 1, 1/2, and 0. 

The monochromatic flux spectrum at the surface of a classical gray atmosphere with net radianl 
flux 77 f will be represented by the symbol l F v (0) and written in the form 


where 


and 


K ( 0 ) 



(147) 


: f 


y 




(148) 


MM 


F 

MM 


(149) 
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Table 4 


The Function f y (t). 



1 

2 

3 

4 

y 5 

6 

8 

10 

12 

14 

0.0 

0.093518 

0.189391 

0.209125 

0.176314 

0.127426 

0.083523 

0.030094 

0.009506 

0.002808 

0.000802 

0.001 

0.093461 

0.189284 

0.209068 

0.176324 

0.127472 

0.083574 

0.030123 

0,009517 

0.002812 

0.000803 

0.01 

0.092952 

0.188334 

0.208561 

0.176412 

0.127874 

0.084024 

0.030381 

0.009615 

0.002843 

0.000812 

0.1 

0.088326 

0.179625 

0.203633 

0.176903 

0.131428 

0.088197 

0.032893 

0.010606 

0.003167 

0.000908 

0.2 

0.083827 

0.171055 

0.198365 

0.176894 

0.134691 

0.092351 

0.035604 

0.011733 

0.003550 

0.001025 

0.3 

0.079822 

0.163354 

0.193302 

0.176474 

0.137402 

0.096096 

0.038238 

0.012884 

0.003956 

0.001153 

0.4 

0.076212 

0.156362 

0.188445 

0.175749 

0.139660 

0.099486 

0.040796 

0.014052 

0.004382 

0.001290 

0.5 

0.072930 

0.149972 

0.183794 

0.174790 

0.141535 

0,102560 

0.043274 

0.015233 

0.004827 

0.001437 

0.6 

0.069927 

0.144108 

0.179345 

0.173652 

0.143081 

0.105347 

0.045670 

0.016422 

0.005289 

0.001593 

0.7 

0.067167 

0.138705 

0.175094 

0.172374 

0.144341 

0.107872 

0.047983 

0.017614 

0.005766 

0.001759 

0.8 

0.064621 

0.133713 

0.171035 

0.170989 

0.145353 

0.110160 

0.050213 

0.018806 

0.006256 

0.001933 

0.9 

0.062263 

0.129088 

0.167159 

0.169524 

0.146148 

0.112230 

0.052358 

0,019993 

0.006759 

0.002116 

1.0 

0.060073 

0.124790 

0.163459 

0.167999 

0.146754 

0.114102 

0.054421 

0.021174 

0.007271 

0.002306 

1.2 

0.056134 

0.117064 

0.156557 

0.164835 

0.147484 

0.117311 

0,058304 

0.023506 

0.008323 

0.002709 

1.4 

0.052691 

0.110322 

0.150264 

0.161598 

0.147695 

0.119905 

0.061873 

0.025787 

0.009400 

0.003138 

1.6 

0.049659 

0.104397 

0.144516 

0.158358 

0.147509 

0.121984 

0.065148 

0.028004 

0.010494 

0.003591 

1.8 

0.046972 

0.099158 

0.139256 

0.155159 

0.147014 

0.123628 

0.068147 

0.030150 

0.011598 

0.004063 

2.0 

0.044577 

0.094498 

0.134432 

0.152030 

0.146284 

0.124907 

0.070892 

0.032221 

0.012706 

0.004553 

3.0 

0.035737 

0.077378 

0.115405 

0.137887 

0.140685 

0.127527 

0.081425 

0.041393 

0.018135 

0.007173 

4.0 

0.030144 

0.066535 

0.102148 

0.126352 

0.134046 

0.126664 

0.088089 

0.048723 

0.023160 

0.009908 
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Figure 1 — The frequency distribution of the net fiux 
of radiation at various levels in an Einstein-gray 
atmosphere. 

where 


B.B.f - 
y 

Values of this function also are listed in Table 5. 


Some numerical values of the function 1 f , 
computed by the same methods used to find 
f (0), are listed in Table 5. These values were 
obtained from two independent determinations 
of the function l z (V), one based on the discrete 
ordinate method (using 96 ordinates), the other 
based on the iterative procedure. 

The monochromatic flux spectrum at the 
surface of the blackbody atmosphere with es- 
caping total flux "tt F will be denoted by the sym- 
bol b b f v (F) and written in the form 

bb -F v (F) = kf^F BB f y . (150) 



Table 5 

Surface Monochromatic Fluxes. 


hv 

f y (0) = 

F „ ( 0 ) 


i f -- 
y 

' X F„ (0)' 

y - kT 

e 

F 

W 

F J 

0 

0.000000 

0.000000 

1 

0.093518 

0.086447 

2 

0.189391 

0.184379 

3 

0.209125 

0.208345 

4 

0.176314 

0.178185 

5 

0.127426 

0.130186 

6 

0.083523 

0.086158 

8 

0.030094 

0.031632 

10 

0.009506 

0.010190 

12 

0.002808 

0.003072 

14 

0.000802 

0.000895 


B.B. f 

y 



0.000000 

0.089618 

0.192817 

0.217847 

0.183875 

0.130577 

0.082653 

0.026458 

0.006991 

0.001635 

0.000351 
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Table 6 lists the maxima of these surface 
monochromatic fluxes and the values of y at 
which these maxima are attained. 

The ratio of the fluxes escaping from the 
gray atmosphere of the second kind to the fluxes 
escaping from the isothermal blackbody atmos- 
phere, expressed as differences in the scale of 
astronomical magnitudes, is 

F„(0) , v 

A n m v - -2.5 log B Bf — • (152) 

If the dimensionless parameter y = hv/kT e is 
used instead of the frequency parameter v, as 
the dependent variable, the magnitude differ- 
erences become 

f y 

A u m y - -2.5 log B ; B 7 f * (153) 

y 

A graph of 4 n m y appears in Figure 2. 

Similarly, the magnitude differences be- 
tween a classical gray atmosphere and an iso- 
thermal blackbody atmosphere with the same 
surface flux are 

T (0) 

A - -2.5 log 57 b- (F y (154) 

or 

*f 

4 I m y = -2.5 Log (155) 

y 

Figure 2 contains a graph of A 1 m . 


Table 6 


Maximum Surface Fluxes „ 


y max 

f 

* y ma x 

1 f 

y max 

B.B. f 

y ma x 

2.77 

0.210698 



2.83 


0.209158 


2.82 



0.218886 



Figure 2 — The monochromatic magnitude differences be- 
tween an Einstein-gray atmosphere and a blackbody at 
the effective temperature T e are given by the graph 
labeled A 11 m ((im^ 11 m = °°). The monochromatic 

y y— r o y 

magnitude differences between a classical gray atmos- 
phere and a blackbody at the effective temperature T e 
are given by the graph labeled A 1 m y . 


The graph of A 1 m y shows that the classical gray atmosphere is brighter than the blackbody for 
y > 5.08 and fainter than the blackbody for y < 5.08. The radiation escaping from the deeper layers 
of an atmosphere is coming from levels at temperatures higher than the surface temperature. 

Since the Planck function increases more rapidly in the ultraviolet than in the infrared, the spectrum 
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escaping from the classical gray atmosphere must be bluer than that of the blackbody, since the 
absorption is "gray." However, because the integrated flux from these two atmospheres is the 
same, the blackbody is brighter at the red end of the spectrum. A detailed discussion of this 
phenomenon has been given by Wildt (References 6 and 7). 

The graph of A 11 m y shows that the gray atmosphere of the second kind is brighter than the 
blackbody for 0 < y < 1.59 and for y > 5.75. It is brighter than the blackbody in the ultraviolet for 
the same reason that the classical gray atmosphere is brighter there. Its excess brightness in 
the far infrared must be due to the high transparency of the gray atmosphere of the second kind 
at small values of v which results from the Rosseland factor 



hv " I 

kT(r ) 


}■ 


This quantity increases with depth and evidently makes for such transparency that the emergent 
flux spectrum is greater than that of the blackbody. Since the integrated flux from the gray at- 
mosphere of the second kind is the same as that from the blackbody, there must be an intermediate 
spectral region where the blackbody is brighter. 


THE LIMB DARKENING IN TOTAL RADIATION: 1(0, /i]/ 1(0, 1] 

The specific intensity of the total radiation escaping from the surface of a gray atmosphere of 
the second kind with effective temperature T e / at an angle 0 (cos" V) to the normal, is 

/•OO 

1(0, M) = I I„ (0, m) dv . 0 < At < 1 . (156) 

•'o 


The limb- darkening ratio, in total radiation, of a gray atmosphere of the second kind is 


1(0, aO 

1(0, 1) 


/•CO 

Jo 


•'0 


l v (0, /*) dv 


I„ (0, l)dv 


0 < M < 1 . 


(157) 


As is well known, the function I„ (0, m) is the Laplace transform of the monochromatic source 
function, i.e., 


l*(0./O 



(158) 
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Since 


T v (T> 


I {' ' “ P [‘ Wf>]} 


dt 


and 


f h. 1 

,r 

j ex P |kT(r )J 

-■] 


l * ( °' M) = { eXP (“ kT(r) " ? J o j 1 ” ex P [' WW]} 


dt I dr . 


(159) 


II v is replaced by (kT e /h) y, hv/kT(r) is replaced by yz(r), and (27r 4 k 4 T e 4 ) / /(i5c 2 h 3 ) is replaced by 

F, 


1(0, M) 


and 


j,, V ( ’ ^ 77 4 /X ^ n; y 3 exp yz(T-) - ^ J {l - exp [-yz(t)]} dt^ dr dy (160) 


1(0, M) 
1 ( 0 , 1 ) 


J | y 3 exp ^yz(T) - i J |l - exp [-yz(t)]j dt^ dr dy 
M J J y 3 exp ^yz(r) - J |l - exp [-yz(t)]} d^ 


(161) 


dr dy 


This limb-darkening ratio was computed for 12 values of/u. The results are listed in Table 7. 

The limb darkening of a classical gray atmosphere will be denoted by I l(0, /i)/ r l(0, l). 
Then 


I l(0 1± ) 
*1(0, 1) 



(0, (J,) dv 


(0, 1) d v 


0 < /x < 1 , 


(162) 
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where 


% (0. /*) 


1_ 


r 

J o 


% ( r r) 


exp 



(163) 


I 



hv 

k J T ( r r) 



(164) 


and 'T ( r r) is the temperature at the optical depth l r in the atmosphere. 


If the variable v is replaced by (kT e /h) y and the function hv/k *T ( r r) is replaced by y I z( I r) the 
limb darkening ratio becomes 


*1(0. M) 
J I(0, 1) 


n oo 

y 3 ^ e *p [y l z (v)] - ij 


exp 


■ tt) d Irdy 


n oo 

y3 | exp [ y Iz ( Ir ^] “ 1 j’ exp^v) 


d I r d x r dy 


(165) 


This ratio was computed at the 12 values of is 
listed in Table 7. 

The limb darkening of the Einstein- gray 
atmosphere is slightly less extreme than the 
limb darkening characteristic of the classical 
gray atmosphere. 

THE FUNCTION 

The angular dependence of the radiation in- 
tensity at the surface of the gray atmosphere of 
the second kind can be written in terms of a di- 
mensionless function &0) , defined by 

1(0, aO = B(0)tf(AO - (166) 

This function is the analogue of the function H(/x) 
introduced by Chandrasekhar into the theory of 
the classical gray atmosphere. 

From Equation 158 it follows that 

I„(0. 0) = B„(0) ; (167) 


Table 7 


Limb Darkening in Total Radiation. 


(j, - cos 0 

1(0, m) 

w 1(0, 1) 

1(0, /z) 
(2) 7 - S —- 
1(0, 1) 

1.00 

1.00000 

1.00000 

0.90 

0.94023 

0.93905 

0.80 

0.88023 

0.87788 

0.70 ! 

0.81993 

0.81641 

0.60 

0.75923 

0.75457 

0.50 

0.69800 

0.69220 

0.40 

0.63599 

0.62909 

0.30 

0.57280 

0.56487 

0.20 

0.50766 

0.49878 

0ol0 

0.43864 

0.42897 

0.05 

0.40082 

0.39087 

0.00 

0.35402 

0.34390 


(1) Einstein-gray atmosphere. 

(2) Classical gray atmosphere. 

Values are based on K. Grossman’s computations of the 
source function of the classical gray atmosphere (using 96 
discrete ordinates). 
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i.e., the specific .intensity of monochromatic 
radiation emerging tangentially from the sur- 
face is equal to the value of the monochromatic 
source function at r = 0. Therefore 

1(0, 0) = B( 0) , (168) 

and 

H(0) n 1 . (169) 

Table 8 lists some values of U(^) obtained 
by computing 1(0, m) from Equation 160 and by 
using the value of B(0) given by Equation 139. 
The computed value of U( 0) differs from the 
theoretical value by 9 x 10" 6 . Table 8 also lists 
comparable values of H(m) taken from Placzek 
(Reference 8). 


Table 8 

The #(aO -Function. 



MO ) 

HO ) 

1.00 

2.824759 

2.90781 

0.90 

2.655919 

2.73059 

0.80 

2.486429 

2.55270 

0.70 

2.316097 

2.37398 

0.60 

2.144650 

2.19413 

0.50 

1.971673 

2.01278 

0.40 

1.796510 

1.82928 

0.30 

1.618032 

1.64252 

0.20 

1.434029 

1.45035 

0.10 

1.239040 

1.24735 

0.05 

1.132214 

1.13658 

0.00 

1.000009 

1.00000 


PHOTON DIFFUSION 


A measure of the progressive reddening that the flux undergoes as it approaches the surface 
is the ratio of the local rate of emission of photons to the local rate of absorption of photons, both 
expressed in units of the dimension [cm -3 sec" *] . This ratio will be denoted as 11 m (r) , while the 
corresponding ratio for the classical gray atmosphere will here be denoted as I m ( l r); the latter 
quantity is the function m(r) of Wildt (Reference 6). 

In the gray atmosphere of the second kind, the number of photons emitted per cubic centimeter 
per second at the level r is 


a(r) 


r f r t-n 

J 0 I 1 ' exp kT(r)J ’ 


B „(T ) 

hv 


dv 


2k 3 T/ r” 

a ( T ) J y 2 e~ yz< ' T%> dy , 


and the corresponding number of photons absorbed is 


a(r) 


1 { 1_exp [ _kT To]} 


J y (T) 
h v 


An expression for (r) is given in Equation 71. In terms of the variable y and the function z(t), 
the number of photons absorbed is 


k 3 T 3 

e 

c 2 h 3 


a(r) 



dy . 
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The function n m (r) is then 


n m(r) 


2 f y 2 e‘ y2 < T >dy 

Jo 





dy 


(170) 


Some values of n m(r) are listed in Table 9. 


Table 9 

Photon Multiplication Factor and Normalized Divergence 
of the Photon Net Flux. 


r 

n,n (T) 

n A(r) 

V 

Mv) 


0.00 

1.16907 

.144617 

0.00 

1.21884 

0.179546 

0.01 

1.16007 

.137984 

0.01 

1.20407 

0.169480 

0.02 

1.15306 

.132741 

0.02 

1.19361 

0.162207 

0.03 

1.14723 

.128338 

0.03 

1.18485 

0.156014 

0.04 

1.14218 

.124478 

0.04 

1.17720 

0.150525 

0.05 

1.13765 

.120996 

0.05 

1.17035 

0.145555 

0.10 

1.11967 

.106877 

0.10 

1.14377 

0.125697 

0.20 

1.09549 

.087168 

0.20 

1.10984 

0.098969 

0.30 

1.07924 

.073424 

0.30 

1.08797 

0.080858 

0.40 

1.06732 

.063075 

0.40 

1.07246 

0.067565 

0.50 

1.05812 

.054932 

0.50 

1.06085 

0.057358 

1.00 

1.03212 

.031122 

1.00 

1.03006 

0.029187 

1.50 

1.02027 

.019868 

1.50 

1.01742 

0.017120 

2.00 

1.01381 

.013620 

2.00 

1.01110 

0.010976 

2.50 

1.00992 

.009821 

2.50 

1.00757 

0.007514 

3.00 

1.00742 

.007362 

3.00 

1.00545 

0.005416 

3.50 

1.00573 

.005693 

3.50 

1.00409 

0.004070 

4.00 

1.00454 

.004517 

4.00 

1.00318 

0.003165 

4.50 

1.00368 

.003662 

4.50 

1.00254 

0.002531 

5.00 

1.00303 

.003023 

5.00 

1.00208 

0.002071 

10.00 

1.00082 

.000816 

10.00 

1.00056 

0.000559 

100.00 

1.00001 

.000013 

100.00 

1.00001 

0.000008 
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From the photon multiplication factor n m(r) it is possible to compute 


- 1 - ^ ' <in) 

This is the divergence of the photon net flux divided by the local rate of emission of photons, a 
dimensionless quantity called by Wildt (Reference 6) the normalized divergence of the photon net 
flux. The corresponding ratio for the classical gray atmosphere will here be denoted as *A (M- 
It is the function A(r)of Wildt (Reference 6). Some values of n A(r) and *A ( I r) are given in 
Table 9. 

The photon flux emerging from the surface of the gray atmosphere of the second kind will be 
denoted by v n D(0) . Clearly, 


n D( 0) = 



dr dv 



dr dy . 


(172) 


The photon flux at the surface of an isothermal blackbody, at the temperature T e throughout, 
can be denoted by 77 b b d(F), where 


BB D( F) = 




dv 


2k 3 T e 3 
c 2 h 3 



(173) 


Then 


n D(0) 

bb D(F) 



e -y2(T) 



[1 



dr dy 


5 ( 3 ) 


(174) 


Computations indicate 


n D(0) 

bb D(F) 


1. 214926 
1.202057 


1.010706 . 


(175) 
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Therefore the radiation from the gray atmosphere of the second kind is in excess of photons by 
1.1 percent when compared with the radiation from the blackbody of the same effective temperature. 

THE RADIATIVE STRESS TENSOR 

Let n n lv (r) , n n 2 v (r) and n n 3 „ (r) denote the principal values of the monochromatic radia- 
tive stress tensor associated with a three-dimensional rectangular coordinate system originating 
at the level t, the first axis of which points in the negative r direction. The integrated radiative 
stress tensor has the components 

II n 1 (r) = [ n n lt „ (r)di/ , n n 2 (r) = 

•'0 


n 2.v ( T ) dv 


and 


n n 3 (r) = f n n 3 v (r)dv . 

•>0 


It is the analogue of the radiative stress tensor of the classical gray atmosphere, which has prin- 
cipal values n j ( r r) , n 2 ( : r), and n 3 ( J r) (Reference 7). The components “llj (t), n n 2 (r) and 
n n 3 (r) are related to the mean radiative stress at the level r, n P(r), by 


“hj (r) + n n 2 (r) + n n 3 (t) = 3 n P(r) . (176) 


The mean radiative stress at the level l r in the classical gray atmosphere is p( : r) (Reference 7). 
The first principal value of the monochromatic stress tensor is 

/•+ 1 

IIn i, v (T) = ^-J I v (r, <&■ ■ (177) 


The corresponding component of the integrated stress tensor is therefore 

r° f +1 

n nj(r) = ^ I* (T, M)M 2 dAidv . (178) 

■'0 •'-l 

Rotational symmetry about the r-axis implies that n n 2 (r) = n n 3 (r). Then from Equation 176 it 
follows that 

n n 2 (r) = n n 3 (r) = \ [3 n P(r) - II n 1 (r)] • 
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(179) 



Since the mean radiative stress at the level r is 


HI 


n P(r) 


2n_ 

3c 


f +1 

I ( T . M) d/idv 


(180) 


n n 2 (r) 


n n 3 (t) 



( r , /x ) d/x dv 


77 

~c 


•'o J- 1 


(181) 


The monochromatic principal values must then be 


"ns.v (T> 


n 3 .v (r) 


77 

~c 


r +1 

I 77 

J (r, — 



(r, /x) m 2 dv 


(182) 


The function “n, v (r) can be computed from the source function B(r). From the equation of 
transfer relevant to the Einstein-gray atmosphere it follows that 


j .co 

E,[t v -T v (r)] 

T v ( T > 


dt 
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r T v (r) 

Bv[tv(t)]E 3 [T v (r) - t„] dt„ 
•'O 


(183) 


Since 


T v ( T ) 


[ { 1_6XP F XT("t)]} 


dt , 


TT 277 

n i.v ( T ) = — 


J |l - exp [- T^Ty B„ ( t )| E 3 1 |l - exp } j 


d£] dt 


IT |l - ex P [” kT(r)] B * ( 4 )J E 3 ^[ j 1 " ex P [" kT(f) | d^ dt . (184) 
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When B v ( t ) is replaced by 



and the resulting expression is written in terms of the variable y - hvAT e and the function 
Z( t ) = [f /B( t )] ? , it becomes 



Therefore 


n n, (r) 



e " y z ( t ) 



dt dy 


+ 



e-y^*) E 3 




(186) 


The second monochromatic principal value, as expressed in Equation 181, is the difference 
of two terms. The first term is 2n/c times the mean intensity J v (t), of Equation 72. The second 
term is 1/2 of the first monochromatic principal value. The difference of these two terms is 
then 


(r) 











(187) 
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Therefore 


n n 2 (r) = 


n n 3 ( T ) 



e -v z (0 



e" yz <£>] d£ 


dt dy 


15 



e-y^t) E 3 



dt dy 


15 



(188) 


The mean radiative stress at the level r is 4tt/3c times the mean intensity J(t), i.e., 


n P( T ) 



(189) 


The difference between the normal component of the stress tensor and either one of the tan- 
gential components can be denoted by n V(r). It is numerically equal to n U 1 (t) - n n 2 (t). Values 
of the ratio 77 V(t)/^ t P(t 7 axe listed in Table 10. For the sake of comparison, values of the ratio 


(i T ) _ n i W " n 2 ( Ir ^ 

P( r r) p( ! r) 


characteristic of the classical gray atmosphere, are also listed in Table 10. 

At all levels within the classical gray atmosphere, the normal component of the radiative 
stress is greater than the tangential components, and t) is a monotonically decreasing 

function of J r which approaches 0 at large values of T r . Within the gray atmosphere of the second 
kind, the normal component of the radiative stress is greater than the tangential component from 
the surface down to r = 1.8. Below this level, each one of the tangential components of the stress 
is larger than the normal component. From the surface to r = 2.8, n V(r)/ 77 P(r ) is a decreasing 
function of r. It then begins to increase, approaching 0 at large values of r. 
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Table 10 


The Difference Between the Normal and Tangential Components of the Radiative Stress 
Tensor Divided by the Mean Radiative Stress. 


T 

n V(T) 

i 



fc - 

K 

H 



T 

P l T 

0.00 

0.337164 

0.00 

0.345794 

0.01 

0.300464 

0.01 

0.307012 

0.03 

0.253873 

0.03 

0.258660 

0.05 

0.222242 

0.05 

0.226000 

0.07 

0.197516 

0.07 

0.200493 

0.09 

0.177093 

0.09 

0.179486 

0.10 

0.168106 

0.10 

0.170270 

0.30 

0.072078 

0.30 

0.073287 

0.50 

0.036717 

0.50 

0.038327 

0.70 

0.019978 

0.70 

0.021949 

0.90 

0.011124 

0.90 

0.013301 

1.00 

0.008285 

1.00 

0,010515 

i.io 

0.006129 

1.10 

0.008385 

1.30 

0.003200 

1.30 

0.005445 

1.50 

0.001438 

1.50 

0.003618 

1.70 

0.000366 

1.70 

0.002448 

1.90 

- 0.000283 

1.90 

0,001682 

2.00 

- 0.000503 

2.00 

0,001401 

2.10 

- 0.000670 

2.10 

0,001171 

2.30 

- 0.000893 

2.30 

0.000824 

2.50 

- 0.001009 

2.50 

0.000586 

2.70 

- 0.001059 

2.70 

0.000420 

2.90 

- 0.001066 

2.90 

0,000303 

3.00 

- 0.001059 

3.00 

0.000258 

3.10 

- 0.001046 

3.10 

0.000220 

3.30 

- 0.001011 

3.30 

0.000161 

3.50 

- 0.000966 

3.50 

0.000118 

3.70 

- 0.000917 

3.70 

0.000087 

3.90 

- 0.000866 

3.90 

0.000065 

4.00 

- 0.000840 

4.00 

0.000056 

5.00 

- 0.000611 

5.00 

0.000013 

10.00 

- 0.000167 

10.00 

0.000000 

50.00 

- 0.000006 

50.00 

0.000000 

100.00 

- 0.000001 

100.00 

0.000000 
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Appendix A 


Additional Numerical Results 


The 

mainder 


following table contains the numerical values of the source function B(t)/f and the re- 
function r (r) at the 602 T-points between r = 0 and r = 100 used in the computation. 
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T 

B(t)/F 

r (x) 


T 

B(t)/F 

r { t ) 

1 

0.0 

C. 443597 

0.624179 

61 

0 .0C03P9 

0.444228 

0.624677 

2 

0.000100 

0.443760 

0.624303 

62 

0.00C398 

0.444243 

0.6 24689 

3 

0.000102 

0.443764 

0.624311 

63 

0.000407 

0.444258 

0.624701 

4 

0.000105" 

0.443768 

C. 624314 

64 

0.000417 

0.444273 

0.624713 

5 

0.000107 

0.443772 

0.624317 

65 

0.000427 

0.444289 

0.624725 

6 

0.000110 

C. 443776 

0.624320 

66 

0.000437 

0.444305 

0.624738 

7 

0.000112 

C.44378C 

C. 624323 

67 

0.000447 

0.444321 

0.624751 

8 

0.000115 

C. 443784 

C. 624327 

68 

0 .000457 

C. 444338 

0.624764 

q 

C.C0C117 

C. 443788 

0.624330 

69 

0.CCC468 

C .444355 

0.624777 

in 

0.000120 

C. 443793 

0.624334 

70 

C. 000479 

0.444373 

0.624791 

11 

0.000123 

0.443797 

0.624337 

71 

0.COO49O 

C. 444391 

0.624805 

12 

0.000126 

C. 443802 

0.624341 

72 

0.000501 

0.444409 

0.624820 

13 

0.000129 

C. 443807 

C. 624345 

73 

0.000513 

0.444428 

0.624834 

14 

0.000132 

C. 443812 

0.624349 

74 

0.000525 

C. 444447 

0.6 24850 

15 

0.000135 

C. 44381 7 

0.624353 

75 

0.000537 

C. 444467 

0.624865 

16 

0 . COO 1 3 8 

C. 443822 

0.624357 

76 

0 • OCO 550 

C. 444487 

0.6 24881 

17 

0.000141 

C. 443827 

0.6243M 

77 

C. 000562 

C. 444507 

0.624897 

18 

0. 000145 

0.443832 

0.624365 

78 

0.000575 

0.444528 

0.6 24913 

19 

0.000143 

C. 443838 

0.624369 

79 

0.000589 

C.44455C 

0.62493C 

20 

0.000151 

0.443843 

0.624374 

80 

0 .000603 

0.444572 

0.624948 

21 

0.000155 

C. 443849 

C. 624373 

31 

0.000617 

C. 444594 

0.624965 

22 

C. 000158 

C. 443855 

0.624383 

8? 

0.000631 

0.444618 

0.6 24983 

23 

C. 000162 

0.443861 

0.624383 

83 

0.000646 

0.444641 

0.625002 

24 

C. 000166 

0.443867 

0.624392 

84 

0.00066 1 

0.444665 

0.625021 

25 

C. 000170 

0.443873 

C. 624397 

8 5 

0.000676 

0 .444690 

0.625040 

26 

C. 000174 

0 . 44 38 8 C 

0.624403 

86 

C.C00692 

0.444715 

0.625060 

27 

0.000173 

C. 443886 

C. 624403 

87 

0.000708 

0.444741 

0.625080 

28 

0.000182 

C. 443853 

C. 624413 

88 

0.000724 

C. 444767 

0.625101 

2Q 

0 • COO l 86 

C.4439CC 

0.624418 

89 

0.CCC741 

0.444794 

0.625122 

30 

0.000191 

C. 443907 

0.624424 

90 

0.000759 

0.444R2 2 

0.625 144 

31 

C. 000195 

C. 443914 

0 .624430 

91 

C.0C0776 

0.444850 

0.625166 

32 

0.000200 

C. 44 3922 

0.624436 

92 

0.000794 

0.444879 

0.625188 

33 

0.000204 

C. 443929 

0.624441 

93 

0.000813 

0.444909 

0.625211 

34 

0.000204 

0.443937 

C. 624448 

94 

0.000832 

0.444939 

0.625235 

35 

0.000214 

C. 443945 

C. 624454 

95 

0.00085 1 

0.44497C 

0.625259 

36 

0.000219 

C. 443953 

C. 624460 

96 

0 .00087 1 

0.445002 

0.625284 

*7 

C. 000224 

C. 443961 

C. 624467 

97 

0.000891 

C. 445034 

0.625309 

3e 

C. 000229 

C. 443969 

0.624473 

98 

0 .000912 

C. 445067 

0.625335 

34 

0.000234 

C. 443978 

0.624480 

99 

0.000933 

0.445101 

0.625361 

40 

0.000240 

C. 443987 

C. 624487 

ICO 

0.000955 

0.445136 

0.625388 

41 

0.000245 

C. 443996 

C. 624494 

101 

0.000977 

0.445171 

0.625416 

42 

0.000251 

C. 444005 

C. 624502 

102 

0.001000 

0.445207 

0.625444 

43 

0.000257 

C. 444015 

0.624509 

103 

0.001023 

C. 445244 

0.625473 

44 

0.000263 

C. 444024 

0.624517 

104 

0.001047 

0.44528? 

0.625503 

45 

0.000269 

C. 444034 

0.624525 

105 

0.C0107? 

0.445321 

0.625533 

46 

0.000275 

C .444044 

0.624533 

106 

0.C01C96 

0.445361 

0.625564 

47 

C.C00282 

0 .444055 

0.624541 

1C7 

0.0C1122 

0.445401 

0.625595 

48 

0 • 00 028 R 

0.444066 

0.624549 

108 

0.001 14R 

0.445443 

0.625627 

49 

0.000295 

C. 444076 

0.624553 

109 

0 • CO 1 17 5 

C. 445485 

0.625660 

50 

0.000302 

C.4440P7 

0.624567 

110 

0.001202 

C. 445528 

0.625694 

51 

0.C00309 

C. 44 40 9 9 

0.624576 

111 

0.001230 

C. 445573 

0.625728 

5? 

C. 000316 

0.444111 

0*624585 

112 

0.C01259 

0.445618 

0.625763 

53 

0. 000324 

0.444122 

0.624594 

113 

0.001288 

0.445665 

0.625799 

54 

0.000331 

C. 444135 

0.624604 

114 

0.001318 

0.445712 

0.625836 

55 

0.000339 

0.444147 

0.624614 

115 

0.001349 

0.445760 

0.625874 

56 

0.000347 

C.44416C 

C. 624624 

116 

0.001380 

0.44581C 

0.625912 

57 

0.000355 

C. 444173 

0. 624634 

117 

0.001413 

0.445861 

0.625951 

58 

0.000363 

C. 444186 

0. 624644 

118 

0 • CO 1 445 

C. 445913 

0.625991 

59 

0.000372 

C. 444200 

0.624655 

119 

0 • CO 1 479 

C. 445966 

0.6 26032 

60 

0.000380 

C. 444214 

0.624666 

120 

0.001514 

0.446020 

0.6 26074 




__ T 

B (t )/F 

r(r) 


T 


rJjJ 

121 

0.001545 

C. 446075 

C. 626117 

181 

0.CC6166 

C. 452588 

0.631227 

122 

C.C01565 

0.446132 

C. 626160 

182 

0.C06310 

0.453153 

0.631372 

123 

C.00 1622 

C. 446190 

0.626205 

183 

0.006457 

0.453403 

0.631519 

124 

0. CO 1 660 

C .446249 

C. 626250 

184 

C.C06607 

0.453616 

0.631669 

125 

C. 00 1698 

C. 446310 

C. 626297 

105 

0.006761 

0.453834 

C. 63182? 

. 126 

0.001738 

C. 446372 

C. 626345 

1 86 

0.006918 

0.454057 

0.631978 

127 

0.00 1778 

C. 446435 

C. 626353 

187 

0.007079 

0.454284 

0.632136 

12P 

0.001820 

C.44650C 

0.626443 

188 

0.007244 

0.454515 

0.632297 

129 

0.001862 

C. 446566 

0.626494 

189 

0.007413 

0.454752 

0.632461 

130 

0 • CO 190 5 

C. 446633 

0.626545 

190 

0 • CO 7506 

C. 454993 

0.632628 

131 

O.C01950 

C. 446703 

C. 626558 

191 

0 • CC 7 762 

0.455239 

0.632797 

132 

C.0019C5 

0.446773 

0.626652 

192 

0.007543 

0.45545C 

0.632970 

133 

C. 002042 

C. 446846 

0.6267C8 

193 

0.C08128 

0.455746 

0.633145 

134 

0.002089 

0.446915 

C. 626764 

194 

0 .008 318 

0.456007 

0.633323 

135 

0.002138 

C. 446955 

C. 626822 

195 

0.008511 

0.456274 

0.633504 

136 

0.002188 

C. 447072 

0.626880 

196 

0.008710 

0.456545 

0.633688 

137 

0.002239 

C. 447151 

0.626941 

197 

0 . CC 891 3 

0.456823 

0.633875 

1 3 P 

0.002291 

C. 447232 

0.627002 

19e 

0*009120 

0.457105 

0.634066 

139 

0.002344 

0.447314 

0.627064 

155 

0.009333 

0.457394 

0.634259 

140 

0.002399 

0.447398 

0*627128 

2C0 

0.009550 

0.457688 

0.634455 

141 

0.002485 

0 • 44 74 8 5 

0.627194 

201 

0.009772 

0.457587 

0.634655 

142 

0.002512 

C. 447573 

C.62726C 

2C2 

0.010000 

0.458253 

0.634857 

143 

0.C0257C 

C. 447662 

C. 627323 

203 

0.010233 

C. 456605 

0.635063 

144 

0. 0026*0 

C. 447754 

0.627398 

204 

0.010471 

C. 458923 

0.6 3 5 27? 

145 

0 • CO 26 92 

C. 447848 

0.627469 

205 

0.010715 

0.459247 

0.635484 

146 

0.002754 

0.447944 

0.627541 

706 

0.010965 

0.459577 

0.635699 

147 

0.00281 « 

0.448042 

0.627615 

207 

0.011220 

0.459914 

0.635518 

148 

0.002R84 

C • 44 PI 43 

C. 627691 

208 

0.011482 

0.460256 

0.636140 

149 

0.002951 

C. 448245 

C. 627768 

209 

0.011749 

0.460606 

0.636365 

150 

0. 00302C 

C.44835C 

0.627846 

210 

0.012023 

0.460965 

0.636594 

1M 

0.003050 

C. 448457 

C. 677926 

1 211 

0.C123O3 

0.461328 

0.636825 

152 

0 . CO 3 1 6 2 

C. 448566 

0.628008 

, 212 

C.C12589 

C. 461699 

0.6 37061 

153 

0.003236 

C. 448678 

0.628052 

! 213 

0.012082 

0.462077 

0.637259 

154 

0.003311 

C. 448792 

C. 628177 

714 

0.013183 

0.462463 

0.6 37541 

155 

C.0033P8 

C .448908 

0.628264 

1 215 

0.013490 

0.462P55 

0.6 37787 

15* 

0.003467 

0.44902B 

C. 628352 

216 

0.013804 

0.463255 

0.6 38036 

157 

0 . 00 35 4 8 

C. 445145 

0.628443 

! 217 

0.014125 

C. 463663 

0.6 38288 

1 5 P 

0 . CO 36 3 1 

C. 449273 

0.628535 

218 

C. 014454 

0.464079 

0.6 38544 

159 

0. CO 3 7 1 5 

0.449400 

0.628629 

, 218 

0 • C 14 79 1 

0.464502 

0.6 38803 

160 

0. C03P02 

0.449530 

C.62e725 

220 

0.C15136 

0.464934 

0 .635066 

161 

0. 003890 

0.449663 

0.628823 

221 

C • Cl 548 8 

0.465374 

0.635332 

162 

0.003081 

0.449758 

C. 628922 

1 222 

C • C 15849 

0.465822 

C. 639602 

163 

0.004074 

C. 445936 

0.629024 

223 

0 .016218 

0.466279 

0.6 39875 

164 

0.004165 

C • 45 CO 7 7 

0.629128 

1 224 

0.016596 

0.466744 

0 .6 4 C 1 52 

165 

0.004266 

C.45C221 

0.629234 

225 

0.016982 

0.467218 

0.640433 

166 

0.004365 

C. 450369 

0.629341 

226 

0.017378 

0.467701 

0.64C717 

167 

0. 004467 

0.450519 

0.679451 

227 

0.017783 

0.468193 

0.641005 

168 

0.004571 

C. 450672 

0.625563 

728 

0 . 0 1 8 1 Q 7 

0.460695 

0.641256 

165 

0.004677 

C .450829 

0.629677 

229 

0 .018621 

0 .4692C6 

0.641551 

170 

0 . 00 4 7 8 6 

C.45C9P9 

0.629754 

230 

0.019055 

0.469726 

0.641P5C 

171 

0.004895 

C. 4511 63 

0.629912 

1 231 

0.019498 

0.470256 

0.642152 

172 

0.005012 

C.45132C 

C.63C033 

232 

0.019953 

C.47C797 

0.642498 

173 

0. CO 5 1 25 

C.45149C 

C. 630156 

1 233 

C. 020417 

0.471347 

0.6 42808 

174 

0.005248 

C. 451664 

0.63028 2 

234 

0.020893 

0.471908 

0.6 4 3121 

175 

0.005370 

0.451842 

0.630405 

235 

0.021380 

0.472479 

0.643439 

176 

0.005495 

0.452023 

0.630539 

236 

0.021878 

0.473061 

0.643760 

177 

0.005623 

0 .452208 

C. 630672 

237 

0.022387 

0.473654 

0.6440P4 

178 

0.005754 

C. 452357 

0.630807 

238 

0 .022905 

0.47425P 

0.644413 

179 

0 • 00 58 8 p 

C.45259C 

0.630944 

239 

C. 023442 

0.474874 

0.644745 

180 

0.006026 

C. 452787 

0.6° 10 84 

240 

0.C235R8 

0.475501 

0.645082 


45 



T 

B(t)/F 

r (x ) 


T 

B(t)/F 

r*x) . 

241 

0,024547 

C • 47614 C 

0.645422 

301 

0.097724 

0.548815 

0.674506 

242 

0.025119 

0.476791 

C. 645766 

30 2 

0.100000 

0.550898 

0.675160 

243 

0.02 5704 

C .477454 

0.646114 

303 

0.102329 

0.553023 

0.675821 

244 

0.0263C3 

C. 478129 

0.646466 

304 

C. 104713 

0.555190 

0.676486 

245 

0.026915 

C. 478818 

0.646821 

305 

0.107152 

0.557401 

0.677158 

246 

0.02754? 

0.479515 

0.647181 

306 

0.109648 

0.559655 

0.677834 

247 

0.028184 

C. 480233 

C. 647545 

307 

0.112202 

0.561955 

0.678516 

248 

0.028840 

C.48C962 

0.647913 

308 

0.114815 

C. 564301 

C. 679204 

249 

0.029512 

C. 481703 

C. 648285 

309 

0.117490 

0.566693 

0.679896 

250 

0.C30200 

C. 482455 

C. 648661 

310 

0.120226 

0.569134 

0.68C593 

251 

0.030903 

C. 483230 

0.649042 

311 

0.123027 

0.571623 

0.681295 

252 

0.031623 

0.484015 

0.649427 

312 

0.125893 

0.574162 

0.682002 

253 

0.032359 

0.484814 

C. 649815 

313 

0.12RB25 

0.576752 

0.682713 

254 

0.033113 

0.48563C 

C.65C209 

314 

0.131826 

0.579393 

0.683429 

255 

0.0338P4 

C • 4 864 6 C 

C.65C606 

315 

0.134896 

0 .582087 

0.684149 

256 

0.034674 

0.487307 

C. 651008 

316 

0.13R038 

0.584835 

0.684874 

257 

0.035481 

C.49817C 

C. 651415 

317 

0.141254 

0.587638 

0.685603 

258 

0.036308 

C. 485045 

C. 651826 

718 

0 . 144544 

0.590497 

0.686335 

259 

0.037154 

C. 489945 

0.652241 

319 

0.147911 

0.593413 

0.687072 

260 

0.038019 

0.450959 

0.652661 

320 

0.151356 

0.596388 

0.687812 

261 

0.C38405 

0.491790 

0.653086 

321 

0.154882 

0.599422 

0.688555 

262 

0.039811 

0.492740 

0.653515 

322 

0.15B489 

C. 602517 

0.689302 

263 

0.040738 

C. 453707 

0.653950 

323 

0.162181 

0.605673 

0.690052 

264 

0.041687 

C. 494694 

0.654389 

324 

0.165959 

0.608893 

0.690806 

265 

0.042658 

C. 455699 

C. 654833 

325 

0.169824 

C.61217R 

0.691562 

266 

0.04365? 

C. 496724 

0.655281 

326 

0.173780 

0.615529 

0.692321 

267 

0.044663 

C. 4^7770 

0.655735 

327 

0.177828 

C. 618947 

0 .6 9 3082 

2 6 P 

0. 045709 

C . 40 8835 

C. 656194 

328 

0.181970 

0.622433 

C. 693846 

269 

0.046774 

0 . 49°9 22 

0.656658 

329 

0.186209 

0.62599C 

0.694612 

270 

0.047863 

C.5C1030 

0 • fc c 7 1 28 

330 

0.190546 

0.629618 

0.695380 

271 

0.048978 

C. 502159 

0.657602 

331 

0.194984 

0.633320 

0.6 96150 

272 

0.050119 

0.503311 

0.658082 

33? 

0.199526 

0.637096 

0.696921 

273 

C. 051286 

C. 504485 

0.658567 

333 

C.2C4174 

0.640948 

0.697694 

274 

0.0624P1 

0.505683 

C. 659053 

334 

C. 208930 

0.644878 

0.698469 

275 

0.05370* 

C. 506905 

C. 6*9554 

335 

C • 2 1 3 79 6 

o.64eeee 

0.699244 

2 76 

0.054954 

C . 50 61 5 C 

C.66Q056 

336 

0.218776 

0.652979 

C.7CC021 

277 

0.056234 

C. 5C542C 

C.66C563 

337 

0.223872 

0.657153 

C.7CC798 

278 

0.057544 

C. 51C716 

0.661076 

33 e 

O.229C07 

0.661412 

0.7C1575 

274 

0.05 8884 

C. 512037 

C. 661595 

339 

0.234423 

C. 665757 

0.702353 

280 

C. 060256 

C. 513385 

0.662120 

340 

0 • 239883 

0.670191 

0.703132 

281 

0. 06 1 66 C 

0.514759 

0.662650 

341 

0.245471 

0.674715 

0.7C3910 

282 

0.063096 

C. 516161 

0.663186 

342 

0.251189 

0.679331 

0.7C4687 

283 

0.064565 

0.517590 

0.663728 

343 

0.257040 

0.68404? 

0.7C5465 

284 

0.066069 

C. 519049 

C. 664276 

344 

0.263027 

0.689849 

0.7C6242 

285 

0.C676CP 

C.52C536 

0.664830 

345 

0.269153 

C. 693754 

0.707017 

2 86 

0. C65183 

C. 522053 

C. 665390 

346 

0.275423 

C. 698760 

0.707792 

287 

0.070795 

C.52360C 

0.665955 

347 

C. 281838 

C. 703869 

0.708565 

288 

0. 072444 

0.525179 

0.666527 

348 

0.288403 

0.709083 

0.709337 

285 

0.074131 

C. 526789 

0.667105 

349 

0.295121 

0.714405 

0.710107 

290 

0.075858 

C. 528431 

C. 667689 

350 

0.301995 

0.719836 

C.71CP75 

291 

0.077625 

C. 530106 

0.668279 

351 

0.309030 

0.725379 

0.711641 

2 Q 2 

0.079433 

C. 531814 

0.668875 

352 

0.316228 

0.731037 

0.712404 

293 

0.081283 

C. 533557 

0.669477 

353 

0.323594 

0.736813 

0.713164 

254 

0.083176 

C. 535335 

C.67C084 

354 

0.321131 

C. 742708 

0.7 13921 

295 

0.C85U4 

C. 537148 

C.67C698 

355 

C. 338844 

0.748725 

C. 714676 

256 

0.087096 

0.538997 

C. 671318 

256 

0.346737 

0.754868 

0.715426 

297 

0. C89125 

C. 54C884 

C. 671944 

357 

C. 354813 

0.761139 

0.716173 

298 

0.09 1201 

C. 542808 

C. 672576 

358 

0.363078 

0.767540 

0.716916 

299 

0.093325 

C. 544771 

0.673213 

259 

0.371535 

0.774076 

0.717654 

300 

0.095455 

C. 546773 

C. 673857 

360 

0.380189 

C. 780748 

0.7 18388 
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T 

b(^)/f 

r ( t ) 


T 

B(t)/F 

r (t) 

361 

0. 389045 

0.78756C 

0.719118 

421 

1.548817 

1.629167 

0.743555 

362 

0.398107 

C. 796514 

0.719841 

422 

1.584893 

1.654781 

0.743523 

363 

0.407380 

C. 8016X5 

C. 720560 

423 

1.62181C 

1.680981 

0.743472 

364 

0.416865 

C. 808866 

0.721273 

424 

i. 659587 

1.7C7782 

0.743406 

365 

0. 426580 

C. 816269 

0.721980 

425 

1.6 58244 

1.735197 

0.743325 

366 

0.436516 

0.823828 

0.722680 

426 

1.737801 

1.763242 

0.743229 

367 

0.446684 

0.831548 

C. 723374 

427 

1.778279 

1.791932 

0.743120 

368 

0.457088 

0.839431 

C. 724061 

428 

1.819701 

1.821282 

0.742956 

365 

0.467735 

C. 847480 

0.724741 

429 

1.862087 

1.851307 

0.742850 

370 

0 • 47 86 3 C 

C. 855701 

0.725414 

430 

1.505461 

1.882025 

0.742707 

371 

0.489775 

C. 864057 

0.726078 

431 

1.549845 

1.913451 

0.742543 

372 

0.501187 

C. 872671 

0.726735 

432 

1.995262 

1 .545603 

0.742366 

373 

0.512861 

0.881428 

0.727383 

433 

2.C41738 

1 .978498 

0 .742176 

3 74 

0. 524807 

C. 890373 

C. 728022 

434 

2.089296 

2.012154 

C .74 1974 

375 

0. 537032 

0.899508 

0.728652 

435 

2.137962 

2.046588 

0.741760 

376 

0.549541 

C.9C8840 

0.729273 

436 

2.187762 

2.CP1820 

0.741535 

377 

0.562341 

C. 518372 

0.729885 

437 

2.238721 

2. 1 17869 

0.741299 

378 

0. 57544C 

C • 52 81 08 

0.730486 

438 

2.290868 

2.154753 

0.741052 

3 79 

0. 588844 

C. 938054 

0.731077 

439 

2.344225 

2.192493 

0.740794 

380 

0. 60256 G 

C. 948214 

0.731658 

440 

2.358833 

2.231110 

0*740527 

381 

0.616555 

C. 958594 

0.732228 

441 

2 .454709 

2.270624 

0.74C251 

38? 

0.630957 

0.969198 

0.732787 

442 

2.511886 

2.311056 

C. 735965 

383 

0.645654 

0.580032 

0.733334 

443 

2.570396 

2.352429 

0.7 3567C 

384 

0.660693 

0.591101 

0.733869 

444 

2.630268 

2.394765 

0.735368 

385 

0.676083 

1.0C2410 

0.734393 

445 

2.691535 

2.438066 

0.739058 

386 

0.691831 

1.013Q66 

0.734904 

446 

2.754229 

?. 482416 

0.7 38740 

387 

0. 707946 

1.025773 

0.735403 

447 

2.810383 

2.527779 

0.7 384 16 

388 

0. 724436 

1.037837 

0.735889 

448 

2.88403 2 

2.574200 

0.738085 

385 

0.741310 

1.050166 

0.736362 

449 

2.951209 

2.621703 

0.797748 

39C 

0.758578 

1.062764 

C. 736821 

450 

3.019952 

2 .6703 14 

0.737406 

3S1 

0.776247 

1.075639 

C. 737267 

451 

3.090295 

2.72006C 

0.737059 

392 

0.794328 

1 .C88796 

C. 7377C0 

45 2 

3.162278 

2.770967 

0.7367C7 

393 

0.812831 

1.102243 

C. 738118 

453 

3.235937 

2.823063 

0.736351 

3^4 

0. 831764 

1.1X5985 

0.73852 2 

454 

3.311311 

2.876375 

0.7 3599? 

395 

0.851138 

1. 130031 

0.73891L 

455 

3.388442 

2.930933 

0.7 35629 

396 

0. 870964 

1 . 144388 

0.739286 

456 

3.467369 

2.986766 

0.735263 

397 

0.891251 

1.159062 

0.739647 

457 

3.548134 

3.043904 

C. 734855 

398 

0.912011 

1.174061 

0.739952 

458 

3.630781 

3.1C2377 

0.734526 

399 

0.933254 

1.185393 

0.740322 

459 

3.715352 

3.162217 

0.7 34154 

400 

0.954553 

1.205066 

0.740636 

460 

3 • 8C 1 894 

3.223457 

0.733782 

401 

C. 977237 

1 .221087 

0 .740936 

461 

3.850451 

3.286129 

0. 793405 

402 

1. COOOOO 

1 .2 37466 

0.741219 

462 

3.581072 

3 .350266 

0.733035 

403 

1.023293 

1.25421C 

0.741487 

463 

4 • C 73 803 

3 .4 159C4 

0.732662 

404 

1.047129 

1.271329 

0.741739 

464 

4-168694 

3.483077 

C. 732289 

405 

1.071519 

1.288831 

C. 741975 

465 

4.265795 

3.551821 

0.731917 

406 

1.096478 

1.306726 

0.742195 

466 

4.365158 

3.622174 

0.731545 

407 

1. 122018 

1 .325022 

0.742359 

467 

4.466636 

3.694172 

0.731176 

4CP 

1 . 146154 

1.343729 

0.742587 

468 

4.570882 

3.767855 

0.730808 

4C9 

1. 174898 

1 .362R5P 

C. 742758 

469 

4.677351 

3.84326 1 

0. 730442 

410 

1 .202 2.64 

1.382418 

C. 742914 

4 7 C 

4.7863C1 

3.920432 

0.73C07P 

411 

1. 230269 

1 .402419 

0.743053 

471 

4.897788 

3.999405 

0.725717 

412 

1. 258925 

1.422872 

0.743176 

472 

5.01187? 

4 .08023? 

0.729359 

413 

1 . 28 8 2 50 

1 .443788 

C. 743282 

473 

5.128614 

4.162947 

C. 729004 

414 

l .31R257 

1.465178 

0.743373 

474 

5.248075 

4.247596 

0.728652 

415 

1. 348963 

1.487054 

0.743447 

475 

5.370318 

4.334226 

0.7 20303 

416 

1. 380384 

1.509426 

0.743506 

476 

5.495409 

4.422981 

0.727959 

417 

1.412538 

1 .532307 

C. 743548 

477 

5.623413 

4.513611 

C. 727618 

418 

1 .445440 

1.555709 

0.743574 

478 

5.754399 

4.60646? 

0.727281 

419 

1 • 47 9 1 OB 

1.579645 

0.743585 

479 

5.888437 

4.701484 

0.726948 

420 

1.513561 

1.604126 

0.743580 

480 

6.C25556 

4.758728 

0.726620 
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T 

B(t)/F 

r (t ) 


T 

B(t)/F 

r (t ) 

481 

6. 165950 

4 .898246 

C. 726256 

541 

24.547089 

17.953264 

0.714667 

482 

6,309573 

5.00C091 

0.725976 

542 

25.118864 

18.359549 

0.714570 

483 

6.456542 

5.104317 

C. 725662 

543 

25 .70395 8 

18.775301 

0.714474 

484 

6,606934 

5.21C978 

C. 725352 

544 

2 6.3C268C 

19.2C074C 

0.7 14380 

485 

6. 760830 

5.32C133 

0.725046 

545 

26.915348 

19.636091 

0.714289 

486 

6.918310 

5.431839 

0.724746 

546 

27.542287 

20.081586 

0.714199 

487 

7, C7 94 5 8 

5.546155 

0.724450 

547 

28.183829 

20.537461 

0.714110 

488 

7.24436C 

5.663142 

C. 724155 

548 

28.840315 

21.003957 

0.714024 

489 

7.413102 

5.782863 

0.723873 

549 

29.51209? 

21.481321 

0.713939 

490 

7.585776 

5.9053 60 

0.723552 

550 

30.199517 

21.969808 

0.713856 

491 

7.762471 

6.03C759 

0.723316 

551 

30.902954 

22.465675 

0.713774 

4 C ? 

7,943282 

6.159067 

0.723045 

552 

31.622777 

22.981188 

0.713694 

493 

8.128305 

6.290372 

0.722778 

553 

32.359366 

23.504619 

0.7 13615 

494 

8. 31 7638 

6.424742 

0.722517 

554 

33.113112 

24.040243 

0.713538 

495 

8.51 13«C 

6.562251 

C. 722260 

c 5 5 

33.884416 

24.588347 

C. 713463 

496 

8.709636 

6.7C2970 

C.7220C8 

556 

34.673685 

25.149220 

C. 713388 

4°7 

8.912509 

6.846974 

0.721761 

557 

35.481339 

25.723159 

C. 713315 

498 

9. 12CI0P 

6.994341 

0.721519 

*5 8 

36.3C7805 

26.310469 

0.713244 

49^ 

9.332543 

7. 145147 

0.721281 

559 

37.153523 

26.911461 

0.713173 

500 

Q. 545926 

7.299473 

0.721048 

560 

38.018940 

27.526455 

0.713104 

501 

9. 772372 

7.457401 

0.720820 

561 

38.9C4514 

28.155775 

0.713037 

502 

10. COOOOO 

7.619014 

0.72C596 

562 

39.810717 

28.799757 

0.7 1297C 

5C3 

10.232930 

7.784399 

C. 720376 

863 

40.738028 

29.4 5P74C 

C.7129C5 

504 

10.471285 

7.953643 

0. 720161 

564 

41.686938 

30.133075 

0.712840 

505 

10.715193 

8.126836 

C. 719950 

565 

42.657952 

2C. 823119 

0.712777 

506 

10, 964782 

8.304069 

0.719744 

566 

43.651583 

31.529239 

0.712715 

507 

11.220185 

8.485437 

0.719541 

56 7 

44.668359 

32.251807 

0.712654 

5CP 

11.481536 

8.671036 

0.719343 

568 

45.7C6819 

32.991208 

0.712554 

509 

11.748976 

8.860964 

0.719149 

569 

46.773514 

33.747834 

0.712535 

510 

12.022644 

9.055323 

0.718959 

570 

47.863009 

34.522086 

0.712477 

511 

12, 302688 

9.254214 

0.718772 

571 

48.977882 

35.314373 

0.712420 

512 

12.585254 

9.457744 

C. 718590 

57? 

50.118723 

36.125118 

0.712364 

513 

12. 882496 

5.666021 

0.718411 

573 

51.286138 

36.954748 

0.712309 

514 

13, 182567 

5.875154 

0.718236 

574 

52.480746 

37.803705 

0.7 12255 

515 

13.489629 

1C .097258 

C. 718065 

575 

53 . 7C318C 

38.672437 

0.712202 

916 

13.803843 

1C .320447 

0.717857 

576 

54.954087 

39.561407 

0.712150 

517 

14.125375 

1C. 548840 

0.717733 

577 

56.234133 

40.471085 

0. 7 1 209 P 

518 

14,454398 

10.782555 

C. 717572 

578 

57.543994 

41.401954 

0.712048 

51° 

14.791084 

11.021726 

0.717414 

579 

58.884366 

42.354507 

0.711998 

520 

15. 135612 

11.266465 

0.717260 

580 

60.255959 

43.329249 

0.711949 

521 

1 5.488166 

11.516918 

0.717109 

581 

61.659500 

44.326697 

0.711901 

52? 

15. 848932 

11.773205 

0.716961 

582 

63.C95734 

45.347380 

0.711854 

523 

16.218101 

12.035467 

0.716816 

583 

64.56542 3 

46.391839 

0.71 1PC7 

524 

16,555865 

12.303842 

C. 716674 

584 

66.C69345 

47.460628 

0.711762 

525 

16.982437 

12.578473 

C. 716535 

585 

67.608298 

48.554313 

C. 711717 

526 

17.378006 

12.855505 

0.716359 

586 

69.183097 

49.673476 

0.711673 

527 

17.78 2794 

13.147087 

C. 716266 

587 

70.794578 

50.8187C8 

0.711630 

528 

16 . 1970C5 

13.441372 

0.716136 

588 

72.443596 

51.990617 

0.7H587 

52° 

1P.62CR71 

13.742516 

0.716008 

589 

74.131024 

53.189824 

0. 7 1 1545 

530 

19.054607 

14.050679 

0.715883 

59C 

75.857758 

54.416967 

0.711504 

531 

19.498446 

14.366024 

C. 715761 

591 

77.624712 

55.672694 

0.711464 

532 

19.552623 

14.688717 

C. 715641 

592 

79.432823 

56.957671 

0.711425 

533 

20.417379 

15.018931 

C. 715524 

593 

81 .283052 

58.272581 

0.711386 

534 

2 0.892961 

15.356841 

C.7154C9 

594 

83.176377 

59.618121 

0.711348 

535 

21.379621 

15.702624 

C. 715296 

595 

85.113804 

6C. 995003 

0.7 11310 

536 

21.877616 

16.056466 

0.715186 

596 

87.096359 

62.403958 

0.711274 

537 

22.387211 

16.41 8553 

0.715078 

597 

89.125094 

63.845733 

0.71123e 

53 8 

22.908677 

16.789078 

0.714972 

598 

91.201084 

65.321092 

0.711202 

539 

23.442288 

17.168236 

0.714868 

599 

93.325430 

66.830818 

0.711168 

540 

23.988329 

17.55623C 

0.714767 

600 

95.499259 

68.375711 

0.711134 





601 

97.723722 

69 . 95659 C 

0.711101 





602 

1CO.CCOOOO 

71.574,29? 

0.711068 
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